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Abstract 

Knowledge of the B s decay fraction of the T(5S) resonance, fs, is important for B s meson studies 
at the T(55) energy. Using data collected by the CLEO III detector at CESR consisting of 0.423 ftr 1 
on the T(5S) resonance, 6.34 fb _1 on the T(4S) and 2.32 fb _1 in the continuum below the T(4S), we 
measure B(T(5S) -► <j>X) = (13.2 ± 0.7+^)% and B(T(AS) -► <f>X) = (7.1 ± 0.1 ± 0.6)%, the ratio of 
the two rates is (1.9±0.1]to'2)- This is the first measurement of the <f> meson yield from the T(5S). Using 
these rates, and a model dependent estimate of B(B S — > 4>X), we measure fs = (27.3 ± 3.2^ 1 gf)%. 
We update our previous, independent measurement of this branching fraction using the inclusive D s 
yields to be (21.8 ± 3.4l|;|)%, due to chan ges in the Df — > (f>ir + branching fraction and a better 
estimate of the number of hadronic events. We also report the total T(5S*) hadronic cross section 
above continuum to be a(e + e~ — > T(56*)) = (0.301 ± 0.002 ± 0.039) nb. This allows us to extract 
the fraction of B mesons as (58.9 ± 10.0 ± 9.2)%, equal to I-/5. Averaging the three methods gives a 
model dependent result of fs = (26^4)%. The results presented in this document are preliminary. 
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I. INTRODUCTION 



The putative T(5S) resonance was discovered at CESR long ago by the CLEO 1] and 
CUSB |2[ collaborations by observing an enhancement in the total e + e~ annihilation cross- 
section into hadrons at a center-of-mass energy of about 40 MeV above the B*B S production 
threshold. Its mass and its production cross-section were measured as (10.865±0.008) GeV/c 2 
and about 0.35 nb, respectively 

The possible final states of the T(5S) resonance decays are: B*B*, B S B*, B S B S , B*B*, 
BB*, BB, BBtt, BB*ti, B*B*tt, BBtttt, where the B and ir mesons can be either neutral or 
charged. Several models involving coupled channel calculations have predicted the cross-section 
and final state composition The Unitarized Quark Model for example, predicts that 

the total bb cross-section at the T(5S) energy is dominated by B*B and B*B S production, 
with B s production accounting for about 1/3 of the total rate. The original ~116 pb~ x of data 
collected at CESR did not reveal if any B s mesons were produced. By measuring the inclusive 
D s production rate using CLEO III data we previously measured the fraction of B s meson 
production at the T(5S), fs, to be (16.0 ±2.6 ±5.8)% of the total bb rate [5j; this measurement 
uses a theoretical estimate of 8(B S -> D S X) = (92 ± 11)% @. 

B s production was confirmed in a second CLEO analysis that fully reconstructed B s 
meson decays. In addition, this analysis showed that the final states are dominated by the 
B*B* S decay channel. These results have been confirmed by the Belle collaboration [7]. A third 
CLEO analysis of the exclusive B reconstruction at the T(5S) showed that the BBX final 
states in T(5S) decays are dominated by B*B* with a considerable contribution from BB* and 
B*B final states 

Knowledge of the B s production rate at the T(5S) resonance, fs, is necessary to compare 
with predictions of theoretical models of 6-hadron production. More importantly, fs is also 
essential for evaluating the possibility of B s studies at the T(5S) using current 5-factories, and 
for future e + e~ Super- B Factories, should they come to fruition, where precision measurements 
of B s decays are an important goal . However, the determination of fs in a model independent 
manner requires several tens of fb _1 • In this paper, we improve our knowledge of the B s 
production at the T(5S) by using inclusive meson yields to make a second model dependent 
measurement of B(T(5S) — > B*B*). 

We choose to examine meson yields because they will be produced much more often in B s 
decays than in B decays. The branching fractions of B s to D s mesons and B to D mesons are 
of order 1, while the rates of B s to D mesons and B to D s mesons are on the order of 1/10. 
We also know that the production rate of mesons is one order of magnitude higher in D s 



decays than in D° and D + decays as measured using recent CLEO-c data 11]. These results 
are shown in Table The rate of D mesons into mesons is only 1%, while the rate of D s 
mesons into mesons is 15.1%. 

Since the B — > 0X branching ratio already has been measured to be (3.5 ± 0.7)% Q, 
we expect a large difference between the yields at the T(5S) and at the T(4S), due to the 

presence of B s , that we will use to measure the size of the B^}B S ^ component at the T(5S). 
The analysis technique used here is similar to the one used in [5j. We also present an update to 
the first measurement of fs that used D s yields and we report on the measurement of the 
total T(5S) hadronic cross-section. When we discuss the T(5S) here, we mean any production 
above what is expected from continuum production of quarks lighter than the b at an e + e~ 
center-of-mass energy of 10.865 GeV. 
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TABLE I: D°, D + and inclusive branching ratios into mesons. The D meson decays used 281 
pb _1 of data at the -0(3770), while the D s decays were measured using 71 ptr 1 at or near 4170 MeV. 



B(D -» 4>X){%) 


D° 


1.0 ±0.1 ±0.1 


D+ 


1.0 ±0.1 ±0.2 


Dt 


15.1 ±2.1 ±1.5 



The CLEO III detector is equipped to measure the momenta and directions of charged 
particles, identify charged hadrons, detect photons, and determine with good precision their 
directions and energies. It has been described in detail previously [3, EH. 

II. MEASUREMENT OF E(T(5S) -> B^TTJ^) USING <j> MESON YIELDS 

A. Data Sample and Signal Selection 

We use 6.34 fb _1 integrated luminosity of data collected on the T(4S f ) resonance peak and 
0.423 fb~ x of data collected on the T(5S) resonance [Ecu = 10.868 GeV). A third data sample 
of 2.32 fb _1 collected in the continuum 30 MeV in center-of-mass energy below the T(4S) is 
used to subtract the four-flavor (it, d, s and c quark) continuum events. 

Hadronic events are selected using criteria based on the number of charged tracks and the 
amount of energy deposited in the electromagnetic calorimeter. To select "spherical" 6-quark 
events we require that the Fox- Wolfram shape parameter R2, be less than 0.25. (ft meson 
candidates are looked for through the reconstruction of a pair of oppositely charged tracks 
identified as kaons. These tracks are required to originate from the main interaction point 
and have a minimum of half of the maximal number of hits. They also must satisfy kaon 
identification criteria that uses information from both the Ring Imaging Cherenkov (RICH) and 
the ionization loss in the drift chamber, dE/dx, of the CLEO III detector. Kaon identification 
has been described in detail previously p. 

B. 4> Meson Yields From T(5S) and T(4S) Decays 

All pairs of oppositely charged kaon candidates were examined for <fi candidates if their 
summed momenta is less than half of the beam energy. Instead of momentum we choose to 
work with the variable x which is the <fi momentum divided by the beam energy, to remove 
differences caused by the the change of energies between continuum data taken just below the 
T(4S), at the T(4S) and at the Y(5S). The K + K~ invariant mass distributions for x < 0.5 are 
shown in Fig. ^ 

1. K + K~ Invariant Mass Spectra and Yields 

A crucial aspect in this analysis is to accurately model the signal and background shapes. 
For the signal we use a Breit-Wigner signal shape convoluted with a Gaussian. The width of 



4 



40000 



> 





C\J 



30000 - 



C 20000 


> 
CD 



O 

=tfe 



10000 - 



I I I I I I I I I I I I I I I I I 

(a) a 








6000 



1630706-016 



5000 - 



4000 - 



- 3000 - 



2000 - 



1000 - 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — T 



(b) 




0.96 1.00 1.04 1.08 1.12 



I i i i I ,, i I i 

0.96 1.00 1.04 1.08 1.12 



> 




w 
c 


> 




1600 - 



1200 - 



800 - 



q 400 
=tfc 



t — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — r 

(c) 




0.96 




l t VW ' l . t *V ' l Vv t l 'r 



■ 11 



1.00 1.04 1.08 1.12 



Invariant Mass (GeV) 



FIG. 1: The invariant mass distributions of the eft candidates with x < 0.5 from: (a) the T(4S) on- 
resonance data (b) the continuum below the T(4S) resonance data, and (c) the T(5S) on-resonance 
data. The solid line is the fit to the signal and background shapes explained in the text. 

the Breit-Wigner function was fixed to the natural width of the <p meson, Y\ = 4.26 MeV [l2T |. 
It is convoluted here with a Gaussian function to allow the integration of the detector resolution 
into the signal function. A second Gaussian is added for an adequate fitting of the tails. The 
form of the signal fitting function, where the dependent variable ( is the invariant mass of our 
K + K~ candidates, is given by 

s(0 = HO + [ C+5a f(y)g((-y)dy. (l) 
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Here, at every physical point £, the Breit-Wigner function 
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at the convolution point y, which is allowed to vary within 5 standard deviations around 
the physical point (. A\ and <Ji are the area and the standard deviation of the corresponding 
functions, and y is the mean value of the Breit-Wigner distribution. h(() is the second Gaussian 
added to fit the tails. 

The background shape is given by a function, chosen to model the threshold, that has the 
following functional form 



where A is the normalization, Co is the threshold, p is the power of a polynomial about the 
turn-on point, and c\ and c 2 are linear and quadratic coefficients in the exponential. Because 
this function has a sharp rise followed by an exponential tail, we are able to accurately describe 
the threshold behavior at low invariant mass close to the kinematic limit. 

We show the invariant mass of the K + K~ candidates in 9 different x intervals (from 0.05 to 
0.50) for all the data samples in Fig. |2l Fig. Eland Fig. 0] We determine the raw yield of 0's 
shown in Fig. |5] (a), (b) and (c) from fits of the invariant mass distributions to the functions 
S(() and -B(C) defined above. All the parameters in S(() in all the fits, except the yields, were 
fixed to the values obtained when fitting the corresponding distributions from the sum of the 
three data samples collected at the three different beam energies in each x interval separately. 
The parameters describing B(Q, on the other hand, are allowed to float except for Co which is 
fixed to twice the K + mass. The yields are listed in the second and third columns of Table HP 
and Table ITTT1 

The systematic errors on these yields are estimated to be ±4%, and are obtained by varying 
the fitting techniques. One variation is to float the fitting parameters at the different beam 
energies, instead of fixing them. We also used different functions including a Breit-Wigner 
signal shape with and without convolution, or a double Guassian shape for the signal, and for 
background either other threshold functions or polynomials. 

2. Continuum Subtraction 

The numbers of hadronic events and <fi candidates from the T(5S) and the T(4S) resonance 
decays are determined by subtracting the scaled four-flavor ( u, d, s and c quarks) continuum 
events from the T(4S) and the T(5S) data. The scale factors, S n s, are determined using the 
same technique described in p] where 



B(() = A • (C - Co) P ■ exp [ Cl (C - Co) + c 2 (C - Co) 2 ] 
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FIG. 2: The K + K~ mass combinations from T(4S) on-resonance data, fitted to the sum of a Breit- 
Wigner signal shape centered at the 4> nominal mass and convoluted with a Gaussian to describe 
the detector resolution, and a second Gaussian distribution for a better parametrization of the tails. 
The background is parameterized by a threshold function (see text). These distributions are in the 
x intervals: (a) 0.05 < x < 0.10, (b) 0.10 < x < 0.15, (c) 0.15 < x < 0.20, (d) 0.20 < x < 0.25, (e) 
0.25 < x < 0.30, (f) 0.30 < x < 0.35, (g) 0.35 < x < 0.40, (h) 0.40 < x < 0.45, (i) 0.45 < x < 0.50. 
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FIG. 3: The K + K~ mass combinations from the continuum below the T(4S), fitted to the sum of a 
Breit-Wigner signal shape centered at the <f> nominal mass and convoluted with a Gaussian to describe 
the detector resolution, and a second Gaussian distribution for a better parametrization of the tails. 
The background is parameterized by a threshold function (see text). These distributions are in the 
x intervals: (a) 0.05 < x < 0.10, (b) 0.10 < x < 0.15, (c) 0.15 < x < 0.20, (d) 0.20 < x < 0.25, (e) 
0.25 < x < 0.30, (f) 0.30 < x < 0.35, (g) 0.35 < x < 0.40, (h) 0.40 < x < 0.45, (i) 0.45 < x < 0.50. 
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FIG. 4: The K + K~ mass combinations from the T(5S), fitted to the sum of a Breit-Wigner signal 
shape centered at the <j> nominal mass and convoluted with a Gaussian to describe the detector reso- 
lution, and a second Gaussian distribution for a better parametrization of the tails. The background 
is parameterized by a threshold function (see text). These distributions are in the x intervals: (a) 
0.05 < x < 0.10, (b) 0.10 < x < 0.15, (c) 0.15 < x < 0.20, (d) 0.20 < x < 0.25, (e) 0.25 < x < 0.30, 
(f) 0.30 < x < 0.35, (g) 0.35 < x < 0.40, (h) 0.40 < x < 0.45, (i) 0.45 < x < 0.50. 
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FIG. 5: Raw (ft yields from: (a) the T(4S) data, (b) the continuum below the T(4S) data, (c) the 
T(5S) data. 

with L n s, L cont , E n s and E cont being the collected luminosities and the center-of-mass energies 
at the T(nS) and at the continuum below the T(4S). We find: 



and 



#45 = 2.713 ±0.001 ±0.027, 
S 5S = (17.18 ± 0.01 ± 0.17) ■ 10~ 2 



(6) 
(7) 
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TABLE II: (ft yields from the T(4S) (N^ 4S)on ), from the continuum below the T(4S) (N l cont ), and 
from the continuum subtracted T(4S) (^Vx(4S))- ^ so listed are the (ft reconstruction efficiencies (e l 



and the partial T(4S I ) 
yields are ±4%. 



4>X branching fractions as a function of x. The systematic errors on the 



Ebeam 



) 



iV T(45)ore 



■^cont 



0.05-0.10 
0.10-0.15 
0.15-0.20 
0.20-0.25 
0.25-0.30 
0.30-0.35 
0.35-0.40 
0.40-0.45 
0.45-0.50 



4938.1 ± 135.2 
22848.9 ± 228.7 
33207.4 ± 269.1 
26020.4 ± 237.6 
15364.1 ± 178.8 

8107.8 ± 128.2 
4263.5 ± 95.9 
2475.3 ± 75.1 

1256.9 ±57.6 



462.4 ±41.0 
2179.0 ±69.0 
3608.7 ±84.1 
3218.7 ±80.7 
2374.2 ± 69.4 
1597.6 ± 57.8 
1113.4 ±48.6 
786.6 ± 42.3 
443.6 ± 30.9 



3683.5 ± 175.1 8.0 

16937.3 ± 295.5 31.9 

23417.2 ± 352.8 47.1 

17288.3 ±323.1 51.3 
8923.0 ± 259.7 49.1 
3773.5 ± 202.5 42.3 
1242.9 ± 162.9 35.8 
341.2 ±137.1 21.6 
53.5 ±101.8 13.4 



1.46 ±0.09 
1.68 ±0.04 
1.57 ±0.03 
1.06 ±0.03 
0.58 ±0.02 
0.28 ±0.02 
0.11 ±0.01 
0.05 ±0.02 
0.01 ±0.02 



TABLE III: (ft yields from the T(5S) {N\^ s)on ), from the continuum below the T(4S) (N l cont are 
the same as in the previous table), and from the continuum subtracted T(5S) continuum subtracted 
(-^-(ss)) - Also listed are the (ft reconstruction efficiencies (e* which are taken to be the same as in the 
previous table), and the partial T(5S I ) — ► (ftX branching fractions as a function of x. The systematic 
errors on the yields are ±4%. 

g^gS N T( 5 S)on Kont ^T(55) g%j g(%) 

0.05-0.10 214.7 ±26.2 462.4 ±41.0 135.3 ±27.1 8.0 2.7 ±0.6 

0.10-0.15 1034.8 ±46.0 2179.0 ± 69.0 660.4 ± 47.5 31.9 3.3 ± 0.2 

0.15-0.20 1432.5 ±53.0 3608.7 ±84.1 812.4 ± 54.9 47.1 2.8 ± 0.2 

0.20-0.25 1151.9 ±47.3 3218.7 ± 80.7 598.9 ± 49.3 51.3 1.9 ±0.2 

0.25-0.30 789.2 ±38.6 2374.2 ± 69.4 390.2 ± 40.4 49.1 1.3 ±0.1 

0.30-0.35 382.3 ±27.9 1113.4 ±48.6 70.0 ±24.1 35.8 0.3 ±0.1 

0.40-0.45 169.7 ± 19.1 786.6 ± 42.3 34.5 ± 20.4 21.6 0.3 ± 0.2 

0.45-0.50 75.0 ± 13.6 443.6 ± 30.9 -1.2 ±14.6 13.4 -0.1 ±0.2 



We estimate the systematic error on these scale factors by obtaining them in a different 
manner. Here we measure the ratio of the number of charged tracks at the different beam 
energies in the 0.6 < x < 0.8 interval. The lower limit on the x interval used here is determined 
by the maximum value tracks from BB events can have, including smearing because of the 
measuring resolution, and the upper limit is chosen to eliminate radiative electromagnetic 
processes. Since the tracks should be produced only from continuum events, we suppress 
beam-gas and beam-wall interactions, photon pair and r pair events using strict cuts on track 
multiplicities, event energies and event shapes. Since particle production may be larger at the 
higher T(5S f ) energy than the continuum below the T(AS), we apply a small multiplicative 
correction of (1.016±0.011)%, as determined by Monte Carlo simulation to the relative track 
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yields. Track counting gives a 1% lower value for S$s and we use this difference as our estimate 
of the systematic error. Note, that the error on the beam energy (0.1%) has a negligible effect. 

These numbers are updated from those reported in (if, due to an increase in Monte Carlo 
statistics, and the use of a more precise calibration of the beam energy. 



3. <p reconstruction efficiency 



The x dependent <ft detection efficiency shown in Fig. [HI is determined by reconstructing and 
fitting candidates in more than 8 million simulated BB generic decays. The low reconstruction 



0.60 I — i — i — r 



0.50 



CD 

"o 

it: 0.40 
LU 

c= 
o 

y 0.30 



CO 
a 
o 

CD 0.20 



0.10 



0.00 



i — r 



i — i — r 



~\ — r 



1 630706-025 

— I — T 



J_ 



0.05 0.15 0.25 

x dpl/E heam ) 



0.35 



0.45 



FIG. 6: 4> reconstruction efficiency from more than 8 million BB Monte Carlo simulated events. The 
errors are statistical only. 

efficiency at small x is due to the fact that as the </> becomes less energetic, it becomes more 
probable that it is formed of a slow kaon (with momentum below 0.2 GeV/c) whose detection 
is very inefficient, since it is likely to be absorbed in the beam-pipe or vertex detector material. 
The efficiency for large x (0.4 < x < 0.5) is somewhat lower than naively expected because of 
the i?2 cut of 0.25. Due to the low efficiency and large backgrounds in the first bin, x < 0.05, 
we do not measure the </> yield in this interval, but will rely on a model of (p production to 
extract a value. 
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4- <p Branching Ratios 



To find the number of hadronic decays produced above the four-flavor continuum, denoted 
as Ny^s)i we multiply the number of events found in the continuum below the T(4S), -/V"x(4S)' 
by the S^s and 6*55 scale factors, and subtract them from the number of hadronic events found 
at each resonance, NyJ nS y. 

N^ s) = N^ 4S) - S 4S * N° f { { s) = (6.42 ± 0.01 ± 0.26) ■ 10 6 (8) 
^T(55) = ^x(55) " S 5S * N°{[ s) = (0.127 ± 0.001 ± 0.016) • 10 6 . (9) 

Using these numbers together with N^ug-) (-^r(5S)) an< ^ e% which are the yield and <fi 
reconstruction efficiency in the i-th x interval respectively, we measure the inclusive partial 
decay rate of the T(nS) resonance into <fi mesons in the i-th x interval as follows 

*vw - **> = J ^ W -^L__ p^a) (10) 

The results are listed in the last column of Table |H] and Table 11111 and shown in Fig. [7| 
Summing the results for 0.05 < x < 0.50 we measure 

B {x>om \T{AS) -> <PX) = (6.82 ± 0.09 ± 0.55)%. (11) 

As the rate in the first interval 0.00< x <0.05 is experimentally not accessible in this analysis, 
we estimate the branching fraction in this interval using our Monte Carlo simulation to model 
4> production. This interval has 3.6% of the total yield, which corresponds to a branching 
fraction of 0.24%. We take the error as equal to the value. The total production rate then is 

B(T(AS) -> (j>X) = (7.06 ± 0.09 ± 0.60)%, (12) 
that when divided by two gives the B meson branching fraction into <fi mesons of 

B(B -> (f>X) = (3.53 ±0.05 ±0.30)% (13) 

which is in a good agreement and more precise than the PDG 0] value of (3.5 ± 0.7)%. For 
the T(5S) we measure for x > 0.05 

B {x>om) (T(5S) -> <j)X) = (12.9 ± 0.7111)%. (14) 

In addition, we find for x > 0.05 

B < T < 8S ^**>-1.9±0.1ijg (15) 



B(T(AS) -> <PX) 

This result is more than 9.6a significant, including both the statistical and the systematic 
errors, thus demonstrating an almost factor of two larger production of <fi mesons at the T(5S) 
than at the T(4S). 
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FIG. 7: The efficiency corrected (f) branching ratio in each x interval from the T(5S) resonance decays 
(the filled circles, in red) and from the T(4S) resonance decays (the open circles, in blue). Errors are 
statistical only, except for the first bin where the T(4S) data point is a theoretical estimate based on 
Monte Carlo, with an error equal to its value. 

5. Determination of The Statistical and Systematic Uncertainties 



Since the measurements presented in the previous section depend on a large number of 
parameters, the corresponding errors on those measurements can have highly correlated con- 
tributions. In this section, we explain the method we used to extract the statistical and the 
systematic errors on our measurements by taking into account the sources of correlations be- 
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tween the different parameters. 

We measure the production rate at the T(4S) and T(5S) using 

i 10 /V* 

B (,>o.05) (T(nS) _ ^ = E 

iV T(nS) ' D l i=2 € 

I ^ Nr(nS)on ~ Sn * N l cont 



I (nb )on " cont (16) 



( N T(nS)on ~ S n * N^ ont ) ■ B X j=2 

The quantities in this equation are defined as: 

• ^r(nS) * s the number of hadronic events at the T(nS) energy after continuum subtraction. 

• A^( n 5) on and N^ ont are respectively the number of hadronic events from the data taken 
at the T(nS) energy and from the continuum data taken at 30 MeV below the T(4S) 
resonance. 

• £>i is the B(4> — > K + K~) branching fraction of 49.1% lil . 

• ^T(nS) i s the number of candidates in the i-th x interval from the T(nS) resonance 
events (i.e. after continuum subtraction). 

• A^( n5 ) on and N % cont are the number of 4> candidates in the z-th x interval from the T(nS)- 
ON resonance data and from the data taken at the continuum below the T(4S) respec- 
tively. 

• e l is the reconstruction efficiency of <fi candidates in the z-th x interval. 

• S n are the continuum subtraction scale factors described in the previous sections. 

The ratio of the <ft production rate from the T(5S) to the <fi production rate from the T(4S) 
is given by 



Aft 

10 JV Y(55) 



H=2 

5(x>0.05)( X ( 4S ) ^ QX) 1 V 10 N - 



JV T(4S) 1 



10 N r ( 4S)on-^* N ion t 



The independent parameters are: S4 and S 5 , e\ Ny, nS ^ m and N % cont . The error on e l includes, as 
shown in Table HV| a 2% error on the tracking efficiency, a 2% error on the particle identification, 
both per track. In addition, it has a contribution of 1% from the error on the fitting method. 
The total systematic error on the <fi reconstruction efficiency is therefore equal to 5.7%, and is 
correlated among all the bins. 

We determine the errors, including all the correlations, by using a Monte Carlo method 
where we generate each independent quantity as a Gaussian distribution using the estimated 
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TABLE IV: Systematic errors on the cf> detection efficiency. 



Systematic error (%) 


Track finding (per track) 


2 


Particle identification (per track) 


2 


Fitting techniques 


1 


Total 


5.7 



mean as the central value and the estimated error as the width. Then we evaluate the relevant 
measured quantities. Applying this method to R gives the spectra shown in Fig. |S1 We do 
this separately for the statistical and systematic errors and then combine them for the total 
uncertainty. This is the probability distribution for this ratio. It is non-Gaussian. 



1630706-020 



80000 



60000 



CO 

Q 



CO 
-Q 
O 



40000 



20000 




R 

FIG. 8: Probability distributions for R smeared by its: (a) statistical error (fitted dots), (b) systematic 
error (filled triangles, in green) and, (c) total error (filled squares, in red). 



We extract the statistical, systematic and total errors on our measurements from the values 
of these distributions corresponding to 68.3% (one standard deviation) of the area under the 
corresponding curves (either statistical, systematic or total) above and below the maximum 
position. We measure 
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R = 1.9±0.ll°| . (18) 


C. Measurement of The T(55) — ► B$ B s Branching fraction 

Eq. (fTHjl and the spectrum in Fig. demonstrate a significant excess of mesons at the T(5S). 

From these results, we can determine the fraction of the Y(5S) that decays into B^B*\ which 
we denote as fs, in a model dependent manner. We assume that the yields at the T(55) 
comes from two sources, B and B s mesons. The equation linking them is 

B(T(55) - <f>X)/2 = f s ■ B(B S - <f>X) + (1 ~ /s) • B(T(AS) - 0X). (19) 

We restrict our consideration to the interval 0.05 < x < 0.50. We obtain fs via the equation 


_ ff(*>o-Q5)(r(5ff) -> ^g) - ^>°- 05 )(T(4g) -> 0X) 
^ S ~ 2£(*>°- 05 )( J B, -f 0X) - i3(^>°- 05 )(T(45) -i- 0X) ' ^ ^ 

The branching fractions B^ >0 ^{T{hS) -> 0X) and £(*>°-°5)(T (45) -> 0X) are given in Eqs. 
ITT1andlT4l It is necessary to estimate B^ x>om \B s — > 0X). We first show that most of the 0's in 
6 decay result from the decay chain B —>■ (D or -D S )X, D ot D s —>■ <pX. The decay rates for B — > 
DX are tabulated by the PDG 0] and reproduced in Table IVl while the value of B — ► -D S X is 
from our previous measurement, updated in the next section. The first line labeled " B product" 
includes measurements of B(B -> D°X) ■ B{D° -> X-7r + ), -6(5 -> L>+X) • £(£>+ -»■ K-n + n + ), 
and B(B D+X) ■ B(D+ — > 07r + ) [11]. We use here the latest D° and D + absolute branching 
ratios ]l7| and the number from S. Stone jll| for Df, to find the B — ► branching fractions 
listed on the next line. CLEO-c has measured preliminary branching ratios for D°, D + and 
mesons into mesons [l]]]. These results are listed in the third line of Table IVl 



TABLE V: B decay inclusive product branching ratios into D°, D + and Df mesons, the branching 
ratios using CLEO-c latest numbers, the inclusive D decay branching ratios into mesons and the 
resulting B — ► (pX product branching ratios. 





D° (%) 


D+ (%) 


Df (%) 


B product 


2.51 ±0.10 


2.16 ±0.11 


0.394 ± 0.007 


B(B -» (D or D S )X) 


66.1 ± 2.8 


22.7 ± 1.5 


11.5 ±1.9 


B((D or D s ) -> <f>X) 


1.0 ±0.1 ±0.1 


1.0 ±0.1 ±0.2 


15.1 ± 2.1 ± 1.5 


B(B -» (£> or D S )X) • or D s ) - 


^0X) 0.66 ±0.10 


0.25 ±0.05 


1.7 ±0.4 



The rates on the last line give the resulting yields from the decay of B into each of the 
individual charm mesons and the subsequent decay of the charm mesons into a 0. The sum 
of these product rates is (2.6±0.4)%. The difference with the measured B — > <pX rate given 
above in Eq. (JT3)l is an additional unaccounted for (0.9±0.5)%. This rate can be accounted 
for from the production of charm baryons or merely by fragmentation processes. The sum of 
B — > DX and D S X branching ratios in the above Table is (100.0 ±3.5)%. We now assume that 
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the number of D plus Df mesons produced in B s decays is the same as in B decays. Using 
our previous estimate of B(B S — > D S X) = (92 ± 11)% |5J, we have (8±11)% of the B s rate into 
charmed mesons that must be accounted for by a mixture of D° and D + , both of which have 
an equal decay rate into 0's. The predicted rates are 



B(B S 
B (x>om \B s 



<j>X) = (14.9 ±2.9)% 
<j>X) = (14.4 ±2.9)%, 



(21) 
(22) 



where we have added in the (0.9±0.5)% from other processes; the rate for 0.05 > x > is taken 
from our Monte Carlo simulation, and amounts to 2.4% of the total yield. 

Solving Eq. (J20)) and using our procedure for finding the errors by generating Gaussian 
distributions for the independent quantities leads to the probability distributions of fs shown 
in Fig. We measure the B^B S ^* ratio to the total bb quark pair production above the 
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FIG. 9: Probability distribution of £(T(55) -» Bf>B\ '), measured using the inclusive yields of 4> 
mesons, smeared by: (a) statistical error (fitted dots), (b) systematic error (filled triangles, in green) 
and, (c) total error (filled squares, in red). 

four-flavor (u, d, s and c quarks) continuum background at the T(5S) energy as 
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B(T(5S) -> B^B { * } ) = (27.3 ± 3.2+^)% . (23) 



D. Estimate of B(T(5S) -> 

We estimate the branching fraction in the 0.00< x <0.05 interval for T(5S) decays by using 
the fraction of <fi meson yield in the first x interval to the <fi yield in the high x intervals. This 
estimate is obtained from a combination of B S B S X and BBX Monte Carlo simulated decays 
at the T(5S) energy. To combine these two types of events, we use the average of our two 
measurements of the f$ fraction, 25%. Here the fraction in the first bin from B decays is 2.6% 
very similar to the 2.4% from B s decays, that makes the error from f$ negligible. Adding this 
estimate to the £(* >a05 )(T(5S) -> <pX) = (12.9 ± 0.7t 2 { 2 4 )% as measured in Eq. JUJ), we find 

B(T(5S) -> 4>X) = (13.2 ± 0.7±H)%. (24) 



III. UPDATED MEASUREMENT OF B(T(55) -> Bf>B s [ *') USING D s MESON 
YIELDS 

Here we update our measurements of the D s yields given in Ref. ||| using our new scale 
factors (Eqs. El and E), and B{D S -> 0vr + ) = (3.5 ± 0.4)% [[J. We have 



B(T(AS) -> D S X) ■ B(D S -> 0tt) = (8.0 ± 0.2+°;^) ■ 10~ 3 (25) 
B(T(5S) -> D S X) • £>(-D s -> 0tt) = (20.3 ± 1.9±|;f) ■ 10~ 3 (26) 
B(T(45) -> L> S X) = (22.9 ± 0.6± 3 |)% , (27) 
B(B^D S X) = (11.5 ±0.3±J;jj)% . (28) 

(29) 

The latter number is in good ag reement with the PDG H3 value of B(B -> D S X) = (10.5 ± 
2.6 ± 2.5)%, based on B{D S -> 0vr+) = (3.6 ± 0.9)%. 
In addition, we find 

B(T(5S) -> D S X) = (58.0 ± 5.5i^)% , (30) 

and 



B(T(5S)^D S X) _ 
B(T(4S)^D s X)- Z - b± ~ - 3 - 
We rewrite Eq. (|19j) for D s mesons rather than for mesons as 



(31) 



B(T(5S) -> D S A) • £(£, -> </»tt)/2 = / s ■ £(5 S -> D S X) ■ B(D S -> 0tt) (32) 

+ • £(T(4,S) -> D S X) • i3(D s -> 0tt) . (33) 



Then using the product of production rates given by Eqs. EB1 and 123 our model dependent 
estimate of B(B S — > D S X) = (92 ± 11)% [5j and the newest estimate of B(D S — > <pir) = 
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(3.5 ± 0.4)% we obtain the B^B S *^ ratio to the total bb quark pair production above the 



(34) 



four-flavor (u, d, s and c quarks) continuum background at the T(5S) energy of 

fs = £(T(5S) -> sWfif) = (21.8 ± 3.4+ 8 ' 5w - 



-4.2) 



The probability distribution of this measurement within the statistical, the systematic and 
the total error is shown in Fig. ITU1 
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FIG. 10: Probability distributions for B(T(55) -> B^B^), obtained by measuring D s yields, and 
smeared by: (a) statistical error (fitted dots), (b) systematic error (filled triangles, in green) and, (c) 
total error (filled squares, in red). 



IV. THE T(55) PRODUCTION CROSS SECTION AND B MESON YIELDS 

We measure the cross-section 

a(e + e- -> T(5S)) = (0.301 ± 0.002 ± 0.039) nb , (35) 

where the first error is statistical and the second is systematic and due to the error in the 
relative luminosity measurement between T(5S) and continuum. 

Measurement of this cross-section allows us to present the previous CLEO exclusive decay 
results for B mesons at the T(5S) 0] in terms of branching fractions. These are given in 
Table IV1I These results are consistent with theoretical expectations [3j, |4j . 
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The inclusive B meson yield is represented as BBX in the Table. This measurement of 
(58.9 ± 10.0 ± 9.2)% is equal to 1-fs under the assumption that there are no non-fefe decays of 
the T(5S). 



TABLE VI: Branching fractions of the possible T(5S) final states either measured directly in this 
paper or estimated using our measurement of the T(5S) total cross section and the cross section and 
upper limits reported in Ref. 



Different components 
of the T(5S) 



WW 



BBX 

B*B* 

BB* + B*B 
BB 

BB { *\ + B^Btt 
BBtttt 



Measured 
cross-section (nb) 



Corresponding 
Branching Fraction 



21.8 ±3.4" 



-8.5 
-4.2 



0.177 ±0.030 ±0.016 
0.131 ±0.025 ±0.014 
0.043 ± 0.016 ± 0.006 

< 0.038 (@ 90% CL) 

< 0.055 (@ 90% CL) 

< 0.024 (@ 90% CL) 



(using D s yields) 
27.3 ± 3.2™ 
(using 4> yields) 

58.9 ±10.0 ±9.2 
43.6 ±8.3 ±7.2 
14.3 ±5.3 ±2.7 

< 13.8 

< 19.7 

< 8.9 



V. CONCLUSIONS 

We measure the momentum spectra of mesons from T(4S) and T(5S) decays, and use these 

data to estimate the ratio f$ of B^b}*^ to the total bb quark pair production at the T(5S) 
energy as (27.3 ± 3.2_ 6 f )%. We also update our previous measurements of the Df rates to 
extract a value of fg = (21.8 ±3. 4142)%. The central value changes mainly because we now use 
B(Df — > (f>7i + ) of (3.5±0.4)%, and the error is lower due to a better determination of the relative 
luminosity. Furthermore, using our cross-section measurement of a^s = (0.301 ±0.002 ±0.039) 
nb for hadron production above 4-flavor continuum at the T(5S) energy, we measure total B 
meson production and extract f$ = (41.1 ± 10.0 ± 9.2)%. Averaging all three methods gives 
fc = (26U)%. 

Our results as summarized in Table IVII1 and at the current level of precision, are consistent 
with the previously published phenomenological predictions 0,0]. 
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TABLE VII: Summary of our results. 



Quantity Measurement 
Results from inclusive cf> measurements 

£(T(4S) -» 4>X) (7.06 ±0.09 ±0.60)% 

B(B 0X) (3.53 ± 0.05 ± 0.30)% 

7 +2.2\qr 



i3(T(55)^0X) (13.2±0.7i 1]( 
B(T(5£) -» <j>X)/B{T{4B) -» <^A) 1.9 ± 0.1±gJ 

>(*)-r(*)\ mq-uqo+W-^ 



B(T(55) -» BY } B^') (27.3 ± 3.2] 



6.1/ 



Results from inclusive -D<, measurements 
£(T(45) -» L> S X) * £(£>, -» 0tt) (8.0 ± 0.2t°;£) • lO" 3 

B(T(55) L> S X) * £(£>, -» 0tt) (20.3 ± 1.9±|;t) • lO" 3 

B(T(4,S) -» L> S X) (22.9 ± 0.6±|:§)% 

B(B->D 8 X) (11.5 ±0.3^|)% 

B(T(5S) -» L> S X) (58.0 ± 5.5± 1 |;|)% 

£(T(55) -» £> S X)/£(T(4,S) -> L> S A) 2.5 ± 0.2±g;| 

g(T(55) -> B^B^) (21.8 ± 3.4±|;|)% 

Results from cross-section and B measurements 

a(e+e~ -> T(55)) (0.301 ± 0.002 ± 0.039) nb 

g(T(55) -» g^gP) (41.1 ±10.0 ±9.2)% 

from combining inclusive 0, D s and I? measurements 
£(T(5S) -> Si*^^) (26±I)% 
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